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We study the time dependence of the grain size distribution N(r, i) during crystallization of 
a d— dimensional solid. A partial differential equation including a source term for nuclei and a 
growth law for grains is solved analytically for any dimension d. We discuss solutions obtained 
for processes described by the Kolmogorov-Avrami-Mehl- Johnson model for random nucleation and 
growth (RNG). Nucleation and growth are set on the same footing, which leads to a time-dependent 
decay of both effective rates. We analyze in detail how model parameters, the dimensionality of 
the crystallization process, and time influence the shape of the distribution. The calculations show 
that the dynamics of the effective nucleation and effective growth rates play an essential role in 
determining the final form of the distribution obtained at full crystallization. We demonstrate that 
for one class of nucleation and growth rates the distribution evolves in time into the logarithmic- 
normal (lognormal) form discussed earlier by Bergmann and Bill [J. Cryst. Growth 310, 3135 
(2008)]. We also obtain an analytical expression for the finite maximal grain size at all times. The 
theory allows for the description of a variety of RNG crystallization processes in thin films and 
bulk materials. Expressions useful for experimental data analysis are presented for the grain size 
distribution and the moments in terms of fundamental and measurable parameters of the model. 

I. INTRODUCTION 

The micro-morphology of solids impacts in an essential way their mechanical, electronic, optical and magnetic prop- 
erties. It is thus an important task to characterize properly the granularity and homogeneity of materials. This 
allows in particular the determination or tailoring of their functionality for the development of new microdevices and 
nanodcviccs. Crystallization and treatment processes define the microstructure of a material. These time-dependent 
processes occurring under adequate thermodynamic conditions are generally spatially inhomogeneous, generating 
grains with different sizes, shape and orientation. 

One of the main physical observable describing the resulting microstructure is the grain size distribution (GSD), 
N(r,t), which is a quantitative determination of how many grains of a certain size are found in the sample at a 
given tinted The vector r is some entity modeling the size and shape of the grains (e.g. diameter of spherical, or 
semi-axes of ellipsoidal grains, number of atoms in the grains, mass of the grains, etc.) and t is the time. The main 
purpose of the present work is to determine and analyze the time-evolution of the grain size distribution obtained 
from the resolution of a partial differential equation (PDE) , describing the crystallization of an amorphous solid in d 
dimensions. The structure of this first-order PDE is motivated by the study of crystallization processes and involves 
two functions T>(r, t) and v(r, t) that describe the source of nuclei and the growth of grains, respectively. 
The crystallization of very different materials leads to distribution laws such as the normal, the Weibull, the Gamma 
(see, e.g., Ref. [If), or most notably the logarithmic-normal (lognormal) distribution^— The shape of the distribution 
depends on the mechanisms involved in the formation of polycrystalline materials. When fitting experimental data 
the choice of the distribution is often not univocal due to the uncertainty of the sampling size considered and the 
precision of the measurements performed^ Nevertheless, it is remarkable that the great variety of crystallization 
processes can be described by only a few distribution laws. This points towards the fact that the detailed knowledge 
of the interactions and mechanisms involved in the crystallization of solids may determine the shape of grains, but is 
not required for the determination of the GSD (see Refs. [l] andliol and references therein.) 

Many theoretical approaches have been proposed to describe the formation of grains during crystallization and they 
can be divided in two main groups: those that describe the process from an analytical point of viewjiSr— mainly in 
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one dimension, and those that are based on numerical first principle calculations ^ 13 ' 16 " — The latter group heavily 
relies on numerical approaches applied to statistical physics. We take the first, complementary approach and derive 
a closed analytical solution to a phcnomenological model that contains the main ingredients of crystallization in d 
dimensions. The PDE and the obtained solution are rather general, but we limit our discussion in this paper to the 
time-evolution of the grain size distribution for crystallization determined by random nucleation and growth (RNG) 
processes. This type of crystallization occurs, for example, in amorphous Si thin films.— i 6 -^ 

The physical picture underlying the present study is motivated by the Kolmogorov-Avrami-Mchl- Johnson (KAMJ) 
thcoryJ^— Starting from a d— dimensional solid in the amorphous phase, we assume that nuclei form randomly and 
homogeneously over the volume of the sample and over time with a constant microscopic rate Iq. Each nucleus 
subsequently grows at constant microscopic rate Vq into a grain until it impinges on other growing grains, at which 
point the growth in the direction perpendicular to the interface stops. We do not include coalescence of grains in our 
model. The central result of the KAMJ theory is the analytical determination of the volume fraction of untransformed 
material Y(t) during crystallization^ (see Sec. Ill C|) . That the KAMJ result is correct within the assumptions made 
is well documented , 22 ' 23 although it was shown to be the large time limit of a more general expression for the transient 
process.— Several interesting extensions and modifications of the theory have also been discussed to account for a 
variety of other crystallization phenomenal 10 ' 23 Here we work within the framework of KAMJ's original model to 
derive and describe the grain size distribution. 

Several groups have studied the grain size distribution during first-order phase transitions taking into account the 
KAMJ result in their considerations^ 2 - - — Most notably, Sekimoto considered a partial differential equation in which 
appears an effective (sometimes also termed actual, transient or average) nucleation rate I(t) = IoY(t) for his study 
of one-dimensional magnetism Similar work followed that was made in the same spirit J^r— This effective time- 
dependence is a direct consequence of the time decay of the fraction of available space for nucleation of the new stable 
phase, derived by KAMJ] 2 ' 10 ' 12 " — On the other hand, the growth rate has been considered constant, v(t) = Vq in all 
these models. 
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FIG. 1: Sketch of the nucleation and growth crystallization process for which the theory is developed. See text for a discussion 
of the figure. Represented are snapshots and distribution N(p,t) (p is the radius of a grain and t the time) at early stages 
(a,d), intermediate stages (b,e), and later stages (c,f) of crystallization. The grain size distribution is calculated exactly for all 
cases in this work. Grey-shaded grains have stopped growing because of impingement, while white grains grow at a rate vq. 



While staying in the spirit of the KAMJ model by considering constant microscopic nucleation and growth rates Iq 
and vq respectively, our approach differs in an essential way with past work on the grain size distribution. Next 
to the effective time-dependent nucleation rate I(t) = IoY{t) first derived by Kolmogorov and Avrami in Refsi 19 ' 20 , 
we also introduce an effective time-dependent growth rate v(t). The physical justification of this time dependence 
not considered previously lies in the fact that the reduced fraction calculated by KAMJ not only reduces the actual 
nucleation rate, but also affects the average growth rate. This is shown schematically in Fig. [TJ The left column 
displays snapshots of the grain distribution found at early, intermediate and later stages of crystallization. Important 
for the present discussion are the shaded grains, which cannot grow further due to impingement. At a given time 
some grains (in white) grow at a constant microscopic rate vq, while others have zero growth rate (shaded grains). 
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In average, the effective growth rate will be less than vq. Comparing Figs. Q] (b) andQJc) we note that the number 
of shaded grains increases with time, implying that the average growth rate decays in time. The right column shows 
schematically the shape of the grain size distribution N(p,t) at these stages of crystallization. The variable p is 
the grain radius. The goal of the paper is to derive analytically and describe, as a function of time and value of 
parameters, this grain size distribution. All the following figures provide exact results for the sketched GSD of Fig. [T] 
The above considerations can be summarized by stating that we place nucleation and growth on the same footing. 
As will be seen in the following sections, the introduction of time dependence is necessary to describe properly the 
crystallization of an amorphous solid. Our model and derivation leads to substantial differences with respect to those 
studied previously. For example, in contrast to previous work our derivation leads to an explicit truncated lognormal- 
type distribution^ In one dimension it is the lognormal distribution. This result is derived, not postulated. The 
time-dependence of the average growth rate is essential to obtain this result. The truncation of the GSD is another 
feature of our model that follows from the effective time dependent growth rate: we obtain a maximal grain size at 
all times. Previous analytical works do not contain such a physical cutoff; they decay to zero only in the limit of the 
infinite grain size. Finally, we obtain an analytical solution of the equations for any dimension d of the crystal, while 
previous analytical considerations were limited to one dimension. This is because we are able to calculate the inverse 
Laplace transform for any dimension d. 

Using the model just described, this paper provides a detailed derivation and generalization to d— dimensions of 
the remarkable result presented in Ref. [24| for three dimensions. We establish classes (subsets) of solutions of our 
partial differential equation one of which leads to a lognormal distribution in the asymptotic limit of large times. The 
proposed determination of lognormal-likc distributions is interesting in itself because contrary to the usual derivation 
it does not rely on a probabilistic argument 1 9 i 17 i 25 i 26 We discuss in depth the conditions under which such distribution 
is found in our model, and how it relates to existing derivations of the lognormal distribution in Sec. IIVE1 We 
emphasize at this point already that the lognormal asymptotic form is not obtained for the general solution of the 
equation we establish, but only for a certain class of functions v(r,t) and T>(r,t). This class is made of functional 
forms known to be relevant for the description of random nucleation and growth crystallization processes and involves 
a time-dependence of the actual nucleation and growth rates. 

To end this introduction, we point out that the results provided in the following sections may be relevant to a variety 
of topics and phenomena. Indeed, the concept of a distribution law is very general and occurs in many fields; for 
example, the fragmentation of solids ; 25 ' 26 gas evaporation^, the distribution of mass in galaxies^ the growth of 
biological tissues^ the distribution of firms as a function of its number of employees^ the production of scientific 
publications as a function of the number of researchers in a group^i etc. All these processes are described in terms 
of a distribution A(r,t), where r must be appropriately defined and is generally a scalar r. In the examples just 
enumerated, the experimental data can actually all be fitted with more or less accuracy by a lognormal distribution 
or a composition thereofi2£ r 30 i 32 Because of the general structure of the partial differential equation discussed here, 
and the fact that the lognormal distribution emerges from fairly general principles, it cannot be excluded that the 
results of the present paper might be applicable to other phenomena in nature, such as those mentioned above. 
This paper is structured as follows. In Sec. [IT] we describe the partial differential equation for the time-dependent grain 
size distribution iV(r, t) and its physical meaning for RNG crystallization models. In this context we introduce the 
functional forms for nucleation and growth rates T>(r, t) and v(r, t), obtained in the context of the KAMJ modelJ^— 
In Sec. Mil we present the general solution of the PDE leading to the central result [Eq. (|2lT) or (|B10a[) ]. We also 
calculate the time-dependence of the maximal grain size [Eq. (|28l) or (|B10c[) ] and provide a simple relation between 
basic measurable parameters of the model [Eq. (|3"Tj) or (jBllI) ]. The general solution is divided into classes according 
to the time dependence of the nucleation and growth rates and the dimensionality of the crystallization process. We 
show that for one specific class, the distribution evolves into a lognormal form in the asymptotic limit t —> oo^ We 
elaborate on the relevance of this result in the following section. In Sec. lIVi we present a detailed analysis of the time- 
and model-parameter dependences of the GSD and discuss our results. We also provide in Sec. IIVEI a discussion 
of our model in the context of other PDEs and derivations of the lognormal distribution. Finally, we conclude with 
Sec. El An example of application of the theory to the analysis of experimental data can be found in Refs. and 
l33l where the case of solid phase crystallization of amorphous silicon thin films has been considered. 
The theoretical calculations are performed with dimensionless quantities. Since one of the purposes of this paper is 
to provide simple closed analytical forms of the time-dependent distribution that can be used to analyze data and 
extract fundamental parameters of a physical system we added two appendixes. Appendix [S] gives a convenient table 
of all parameters, variables and relations between constants defined in the paper, and Appendix Sec. |B] summarizes 
the main results written in quantities with physical dimensions. 
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II. MODEL 



The theory introduced in this section is potentially applicable to a variety of phenomena in nature. Therefore, we 
first describe in part fSec. Ill A)) the differential equation for the distribution N(r,t) in general terms. We then write 
in Sec. Ill Bl the equation specifically for RNG crystallization. 



A. Differential equation for the distribution N(r,t) 

We are interested in phenomena describable by a time-dependent distribution N(r,i) of a certain entity defined in 
terms of a d— dimensional vector r. In the case of crystallization processes, the entity is the crystalline grain and the 
vector r describes the shape and size of a grain. For example, if the entities are grains of ellipsoidal shape, they may 
be described by the semi- axes rj (j = 1, . . . , d) and r = (ri, . . . , r^). Another possible choice of r is to describe the 
grain in terms of the number of constituting atoms or the mass M, in which case r = M is a scalar (as discussed in the 
next section, the latter choice leads to a non-linear differential equation, which cannot be solved analytically, while 
the first choice leads to a set of linear differential equations that can be solved analytically). These two examples show 
that, in general, the dimension d of the vector r may be different from the spatial dimension. In the considerations of 
the present paper the two dimensions will turn out to be identical. Note that this definition of N(r,t) presupposes 
that the distribution is spatially homogeneous and does, therefore, not depend on the position vector x in the sample. 
We assume that the dynamics of the distribution is affected by two processes. One is a source-and-sink term T>(r, t) 
that describes the creation and/or annihilation of the entities. The other describes the dynamic growth and/or 
shrinking of the entities and is defined by the vectorial quantity v(r, t). The analytical calculations of this paper can 
be extended to PDEs that include further terms to account, for example, for the possible coalescence of entities, as 
long as the equation remains first order and linear. Beyond that, numerical techniques are likely the path to follow, 
except in one dimension, as shown in Refs. [j^-[l4l 

The contribution of the source-and-sink term to the time-evolution of the grain size distribution N(r,t) is given by 

^=X>(M). (1) 

In the following we consider the separation of variables T>(r, t) = I(t)D{r). The function D(r) may often be modeled 
by a Gaussian centered on some characteristic quantity r c (for example, the radius of a spherical nucleus p c ), which 
may be replaced by a Dirac delta distribution in the limit of zero variance. The time-dependent part of the nucleation 
rate, I(t), will be discussed in more details in Sec. Ill Cl below. 

The growth of a fixed number of entities is described in terms of a general continuity equation in the space of the 
vector r 

^M+V r -[JV(r,i)v(r,t)]=0. (2) 
The growth rate of the grains, v(r, t), will also be discussed in more details in Sec. Ill CI 

Assuming that these are the two only contributions to the rate of change in the distribution N(r, t), we are led to the 
following partial differential equation 

+ V r • [A(r, t) v(r, t)} = 2?(r, t) . (3) 

We emphasize the point made above that r is not in general related to a spatial coordinate, but characterizes the 
entities that are created and grow in time. 

The PDE introduced here on physical grounds is well known and displays a large variety of more or less complicated 
solutions, depending on the ingredients introduced to account for a particular phenomenon. The study of this equation 
leads to an interesting and unexpected outcome when applied to crystallization processes, as shown in later sections 
(see also Ref. HJ.) A formal general solution can be found for this first-order linear partial differential equation and 
specific solutions can be obtained by quadrature for well-behaved source-and-sink terms 2?(r, t) and growth rates 
v(r, t)i24 We discuss in the next sections closed analytical solutions for growth rates and source terms applicable 
to random nucleation and growth processes during crystallization of amorphous solids. Note that altered forms 
of the above equation, some pertinent to other physical situations, have been discussed in the literature (see, e.g., 
Refs. [Tllfl2ll35M36l Il3l. and IbH and references therein.) The choice of v(r, t) and T>(r, t), the formalism, and the resulting 
analytical solutions distinguish the present work from others. 
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B. Differential equation for RNG processes 



To apply the differential equation of the previous section to RNG processes, we need to specify the rates T>(r, t) 
and v(r, t). Transmission electron microscope 4 ^— ' 29 ' 37 or atomic force microscope^ measurements performed during 
crystallization of many different materials reveal grains with a great variety of shapes and sizes. These grains are 
generally characterized by a scalar quantity, such as the number of atoms, the volume or the mass of the grain. 
In the present paper, we take an alternative approach, and describe the grains as ellipsoids with semi-axes rj with 
j = l,...,d. For thin films d = 2, whereas for bulk solids d — 3. The components of the vector r defined in the 
previous section are now given by the semi-axes of the ellipsoid r = (?"!,•• • , r d ) with n > • • • > r d > 0. In this 
representation the dimensionality d of the vector r coincides with the spatial dimension of the nucleation and growth 
process. 

In general, the experimental representation of the grain size distribution is done in terms of the average radius of 
a d— dimensional sphere, the volume of which equals the measured volume of the grain. For this reason, we later 
consider the special case of spherical grains to allow comparison of the theory with existing experimental data i 24 ' 33 
This reduction to spheres in d dimensions has the additional advantage of expressing the grain size distribution in 
terms of two variables only (the radius of the grain p and time t) and allows a better understanding of the behavior of 
the distribution as a function of nucleation and growth rates without the additional difficulty related to the anisotropic 
growth and orientation of the grains. 

In our RNG model, the source T>(r,t) describes the formation of nuclei (we do not consider the dissolution and 
coarsening of grains) that is determined by microscopic interactions and the thermodynamical conditions, under 
which crystallization of a specific material occurs. We assume that nuclei are formed with a critical volume Q Ctd at a 
rate I(t). The source term thus takes the form 

V(r) = I(t) D(r) = I(t) 5 (Q d - fi c>d ) , (4) 

where 

n d/2 

&d = w d r x . . .r d , with uj d = —j2 -y, (5) 

gives the volume of a d— dimensional ellipsoid, and T is the gamma function. fl Ctd is defined as in Eq. ([5]) and is the 
critical volume of a nucleus, that is, the volume of the smallest grain that can be found in the sample. This simplified 
expression for the nucleation term which states that all nuclei are formed with the same volume fi Cj d, is consistent 
with the fact that we do not discuss the details of the nucleation process. It is a good first approximation, except for 
describing early stag es of crystallization.— This will become apparent in Sec. IIVI and the application of the theory 
presented in Ref. l33l 

Once a nucleus is formed, crystallization leads to its growth into a grain. One simplification relevant to some RNG 
processes is the case where the variation in the growth rate v(r, t) only weakly depends on r. In this case, Eq. J3]) 
can be approximated by 

^^l + v (t)-V r N(v,t)=V(r,t). (6) 

In fact, this equation applies to the case of the solid phase crystallization processes considered in Refs. andl33l In 
this case, the growth rate v(r, i) is independent of r. This non-trivial simplification is one reason for choosing the 
present formalism in the vector space of r rather than writing the grain size distribution N(y, t) in terms of a scalar y 
such as the volume or the number of atoms of a grain. Indeed, in processes such as solid-phase crystallization^— & 24 i 33 
atoms are present in the immediate vicinity of nuclei and grains at all time. Thus, their contribution to the growth of 
grains does not require diffusion processes and the rate v does not depend on the size of the grains* 1 ^ - — It is essential 
to realize that writing the equation for random nucleation and growth, in terms of a scalar, such as the number of 
atoms in the grain, would result in a non-linear dependence of the growth rate on that scalar, which is generally 
untractable analytically. The present formalism allows to circumvent this problem at the cost of a multidimensional 
PDE, which turns out to be solvable analytically. 

As mentioned above, application s 24 ' 33 assume grains of spherical shape and the grain size distribution depends only 
on the radius p of the d— dimensional sphere. In the space of semi-axes of the ellipsoid, the vector for a sphere is 
r = (fx, ... , r d ) = (p,...,p), and the magnitude of the vector appearing in Eqs. is r = \fdp. Introducing polar 

coordinates, where p is the radial coordinate, and using the expression for the source term ((5]), we obtain 



(7) 
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where A c ^ = uiddp^T 1 is the surface of the hyperspherical nucleus of radius p c in d— dimension. 

We derive in the next section the exact solution of the above equation for rates relevant to RNG processes, with the 

boundary conditions 

N(p<p c ,t) = 0, N(p^ oo,i) = 0, for t>* , (8a) 

and the initial condition 

JV(p,*)=0, for * < t , (8b) 

where t is the incubation time. We emphasize that contrary to many approaches, the present formalism does not 
assume prior knowledge of an initial distribution; there are no grains in the samples until the incubation time is 
reached. 



C. Random nucleation and growth rates 

We define in this section the effective nucleation and growth rates I(t) and v(t) that appear in Eq. ([3]). The thermo- 
dynamics of nucleation and early stages of crystallization is involved and a subject of its own . 10 ' 39 We assume that 
the conditions necessary for the creation of nuclei and their subsequent growth into grains are fulfilled. The creation 
of nuclei is described in terms of the existence of a critical average volume C and an effective rate I(t) introduced in 
Eq. ([4|. We are mainly interested in the time dependence of the grain size distribution, resulting from RNG processes 
and consider situations where no coarsening occurs.^*^ Full crystallization has been completed once each atom of the 
sample is assigned to a grain. This is a realistic scenario, for example, in solid-phase crystallization of an amorphous 
sample; 24 ' 33 and is the one considered in the KAMJ theory for RNG processes ^ 2 1 ' 24 

The model assumes homogeneous nucleation; nuclei appear randomly in the sample (no pre-existing nuclei or nu- 
cleation centers) with a constant microscopic rate Jo, irrespective of the presence of already transformed material. 
To satisfy this condition, the concept of phantom nuclei was introduced in Ref. I20I . Avrami obtained an expres- 
sion relating the effective to the extended volume fraction of transformed material that corrected for the presence of 
these phantom nuclei and lead him to derive an expression for the fraction of material available at time t for further 
nucleation and growth^— 

to is the incubation time, t c is the critical crystallization time, and the integer n determines the time dependence of 
the crystallization process. 0(t) is the Heavisidc function. 

Later work on grain formation in various systems to which the KAMJ applies (e.g. Refs. [Iclfl2l - [l8l and H^) has 
recognized that though the microscopic nucleation rate Jo is constant, the actual (effective) nucleation rate relevant 
for the grain size distribution decays in time as I(t) = IoY n (t) because of the reduced fraction of available space 
for nucleation. This effective rate has been confirmed, refined, and extended to include a variety of nucleation 
phenomena. For example, the expression for I(t) has been generalized to account more precisely for the early stages 
of nucleationJ^ The concept of "phantom nuclei" introduced by Avrami 2 ^ and the importance of which has been 
thoroughly discussed 2 ^ has been linked to the concept of spatial correlations of nuclcation i 12 i 22 Further extensions of 
the model have been proposed to include simultaneous nucleation, non-random processes, coalescence, or to take into 
account the symmetry of the crystal structure in the transformed phase. 12 i 14 ' 22 i 23 

In the present paper we account for the insight gained in earlier work, but extend the analysis in an essential way 
by also considering the effective change in the growth rate. This aspect of the growth of grains during crystallization 
has not been addressed in any previous work. We start with the random and homogeneous nucleation of grains in 
the sample. As the crystallization process develops, each nucleus grows into a grain at a rate vq. The growth of each 
grain eventually comes to a halt, cither because it is inhibited by the presence of another grain boundary (e.g., by 
impingement) or because crystallization has been completed (no further free atoms are available). Thus, in the same 
spirit as for the analysis of the effective nucleation rate, the growth rate actually decreases with time (see Fig.Q]). The 
proper account of this effective time-dependent growth rate v(t) turns out to modify, in important ways, the grain 
size distribution and leads to the lognormal-like form of the distributions observed in experiment iSi^ 
Based on Eq. (|9|), we thus develop further the theory of RNG by introducing different rates I(t) and v(t) for nucleation 
and growth, respectively. We postulate that both rates take similar functional form, except, for different critical times 
t c i and t cv and different power laws in the exponential. As a result, we introduce the following expressions for 
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nuclcation and growth rates: 

In{t) 
V m (t) 

While the effective nucleation rate [Eq. (|10a|) ] has been used in previous work-i2~— to determine the grain size distri- 
bution, the introduction of the effective growth rate [Eq. (jlObp ] is new. The analytical form of v m (i) and the choice 
m = (exponential time decay) will be justified a posteriori by analyzing the analytical results and comparing the 
theory with experimental data<22 It is important to realize that the KAMJ model relies on the fundamental assump- 
tion of constant and homogeneous microscopic nucleation and growth rates, denoted Iq and Vq. By introducing the 
above effective rates [Eqs. (|10p j. we thus presuppose that the physical picture underlying the present theory is the 
same as that of the KAMJ model, but we account for the effective time decay of the nucleation and growth rates. 
The nucleation and growth rates contain five parameters: Io, vq, to, t c i and t cv . However, if one shifts the time 
variable by to, it turns out that only two parameters actually affect the normalized distribution, namely, vo and the 
ratio t cv /t c i. In addition, the integers n and m determining the power law of the time decay have to be specified. For 
m = (or n = 0), the decay is exponential, and for m = 1 it is Gaussian. For m > 1 it is super Gaussian. For m = — 1, 
the growth rate is constant, U-i(t) = evo and the PDE for d = 1 reduces to that of Refs. [ID and [Til in the absence 
of coalescence. The nucleation rate I(t) is essentially given by the fraction of material available for crystallization 
and thus n = d, the dimensionality of the grain^ Consequently, our t c j is the critical time t c in the conventional 
KAMJ model. On the other hand, the growth rate v(t) is determined by another mechanism, namely, the inhibition 
of the grain growth by their neighbors. We therefore expect the exponent m to differ from n. Of importance for the 
following are the values m — 0, 1 and n = 1,2,3. The case m = — 1 and n = 1 has been studied analyticall y 12 i 14 
taking into account coalescence, which is absent in the present model. Note finally that m must not necessary be an 
integer, but for simplicity we consider only this case here. 



IoY*(t) = Iocxp< 
v oYZ(t) = i> exp 
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III. SOLUTION FOR THE EXTENDED KOLMOGOROV-AVRAMI-MEHL- JOHNSON MODEL 

We calculate the grain size distribution N(r, t) for the random nucleation and growth crystallization model described 
in the previous section by solving Eq. ([3]) with Eqs. (|10[) . This can be achieved formally by various methods as, for 
example, the Laplace transform or the methods of characteristics^ An explicit expression for N(r, t) is then obtained 
by quadrature once v(r, t) and T>(r,t) are defined. For the nucleation and growth rates defined in Eqs. (fTU)) . it is 
more natural and instructive to use the Laplace transform. Remarkably, we can obtain the inverse Laplace transform 
in any dimension within the present model. The details of the derivation of N(r, t) are presented in Appendix [Cj We 
also determine the time-dependent maximal size of the grains that can be found in the sample, and the moments of 
the distribution. 



A. Partial differential equation in dimensionless quantities 



It is useful to introduce the following dimensionless variables: 



yjt cv t c j 

and constants 

Zb = iciIo/C V = ^. (lib) 

It is important to bear in mind that a simplified notation is used in Eqs. (|lip and the following to avoid overloading 
expressions with indexes. The dimensionless variable 7 actually depends on m through pi^. This choice of dimension- 
less radius is made because, in general, p c <ti p 1 ^ and the latter is easily measured experimentally while the former 
is not. Since we consider solutions obtained for different growth laws (distinguished by m) separately, no confusion 
can arise from this simplified notation. In Sec. IIV1 we analyze the dependence of the general grain size distribution 
on the parameters of the model, that is, for various values of the constants defined in Eq. (jTTbJ) . 
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To solve the partial differential equation ([7]) it is useful to write it in dimcnsionless variables (jlll) and to introduce 
the auxiliary dimcnsionless function N("f, r) = (p^) 2 A 00 ^ d -f d ~ 1 N(-f,T), where A^^ = Wdd{p'^ l ) d ~ 1 is the surface 
of the largest hyperspherical grain of radius p^ found in a d— dimensional sample at full crystallization. Then, the 
differential equation. ([3]) or Eq. ([7]) takes the form 

+ — — = tr In{T) 6(1 ~ 7c) ' (12) 

with 7 C = Pc/Pmi an d the dimensionless nucleation and growth rates 

J„(r) = 2oFi(r) =2 exp (r - r )" +1 } 6 (r - r ) , (13a) 

«™(t) = V r;;(r) = V cx P j - + j 6 (r - r ) , (13b) 

It is important to remember that the grain size distribution is given by iV(7, r), not by N(~/,t). The latter is only 
an auxiliary function to bring the partial differential equation in a form that can be solved analytically. 



B. Solution of the partial differential equation (112[) 



We derive in Appendix [C] the solution of Eq. (|6]), written in dimensionless quantities as Eq. (fT"2)) , for the KAMJ model 
[Eq. ©], and for effective nucleation and growth rates [Eqs. (jT3")) ]. We obtain the following grain size distribution for 
any dimension d 



N(j,t) 



N(rf,r) 



(p£) 2 4»,d7 
— Z^ ex P 



d-l 



(14) 



7' 



m+1 



- r ) 



n+1 



0(r - cr m ^) 9 {a m<i - r ) 



with the constant prefactor Cd defined by 



«0/ Axj,d' 



(15) 



with Aoo : d defined above Eq. (TT^T) . 

The general solution is directly proportional to the ratio Iq/vq and inversely proportional to r y d ~ 1 . The coefficient 
Cd ~ Iq/vq cancels out when considering normalized distributions. Such is the case in the numerical calculations of 
the next sections and often also in the description of experimental data^ The functions ct„ v ;(7, r) (i = 1,2,...) are 
solutions of (see appendix [C]) 



with u m defined by 



7 = 7c + U m (<7 m .i,T), 



1 f b 

u m (a,b) = —l v m (r)dT. 

tr J a 



(16) 



(17) 



The function <j m ^ only depends on the growth rate, not on the nucleation rate. For the KAMJ growth rate (|13b[) . 
the latter equation has a. single solution; the sum over / in Eq. (1141) and the index i in a,,,,; can therefore be dropped. 
Equation (JT5J) can be rewritten in the form as 



r 



1 



TO. + 1 ' 



To 



tr 



m+1 



(TO + 1) 



7 - 7c 
V 



1 



TO + 1 ' V t. 



m+1 



for m > 0, where r[a, x] is the upper incomplete gamma function. This equation yields for m = 

-l 

00(77 t) = t + t r hi 



7 ~7c 



V 



+ e -(-r-T )/t r 



(18) 



(19a) 
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and for m = 1 



o-i (7, t) = t + t r erf 1 



erf 



T ~ Tp 



2 7- 7c 



V 



(19b) 



where erf is the error function. For higher values of m, Eq. (|16[) is implicit and solving it involves special functions. 
To bring to the fore the physics contained in our general solution for the KAMJ model, we transform the time- 
dependent product of Hcavisidc functions appearing in Eq. (|14j) into a size-dependent difference of these functions 
(see Appendix [C]) 



Q[t 
lim QW 



t (7,r)] e[<7 ro (7,r)-ro] = 9 (7 - 7c ) - 6 [7 - 7^ 
t (7,r)] e[<7 m (7,T)-ro] = 9 (7 - 7c ) - 9 (7 - 1) 



(20a) 
(20b) 



where 



v- 

17) 



; (r) is determined below. This difference expresses the fact that at all times, only grains with radius 



between the minimum p, 
Eq. (TH)) now reads as 



(7 = 7c) and the maximal size p™ ax (i) (7 = 7™ ax ( T )) can be found in the sample. Thus, 



AT(7,r) 



C d 



m+l 

a m T \ , n +l / „_ \n+l 



:{e(7-7c)-e[7- 7 r x M]}- 



(21) 



This is the central result of the theory, the general solution of Eq. (|12l) for spherical grains in d dimensions using 
random nucleation and growth related rates that contain the fraction of material available for further crystallization 
of the Kolmogorov-Avrami-Mehl-Johson model [Eqs. (|T5|) ]. The function a m (^ : T) appearing in Eq. (pH]) is a solution 
of Eq. ([16)) and is given by Eqs. (p~9|) . The maximal grain size 7™ ax is obtained in Eq. (|28[) below. 
This result has several general properties. First, the solutions can be divided into classes characterized by the 
dimensionality of the crystallization process d, and the time-decay of the nucleation and growth rates, which are 
specified by the values of n and m. Thus, each triplet of non-negative integers (d,n,m) defines another class of 
the general solution. In the RNG process discussed in the present paper, n is identified with the dimension d and, 
therefore, the classes are defined by the doublet (n = d, to). Second, the time dependence of Eq. (|2"Tj) appears through 
cr m (r) [Eq. (TTS])] and 7^ ax (r) [Eq. ([2"5])]. which are non-trivial functions of r. Third, the nucleation rate Iq appears 
explicitly only in the prefactor Cd, whereas the growth rate vp is present in a m and 7™ ax as well. This is also the 
case for m = —1 (constant growth rate) and agrees with Ref. [l2|. On the other hand, both critical times t cv and t c j 
appear in the exponential through the ratio t r . Fourth, the derivation of the grain size distribution provides cutoffs 
at the radii of the nucleus and the largest grain found in the sample at time t. This will be discussed further in the 
next section. Finally, we emphasize that the explicit analytical solutions were derived here for specific source and 
growth terms v(r, t) and T>(r,t), namely, Eqs. ([TO]) , which are consistent with the KAMJ model. This is important 
to remember for the discussion of the next sections, especially the result obtained in the limit t — > 00. 
Explicit expressions of the GSD can be written for specific classes of solutions (d,n,m), in particular, for the cases 
of interest for the description of experimental data, namely, m = 0, 1 and n = d = 1,2,3. For m = but arbitrary 
(<i, n), Eq. (f2"Tj) reads as 



iV( 7 ,r) 



Cd 



exp{(-l)« [^lnao]" +1 } 



Of) 



x {0 (7 - 7c) - 9 [7 • 



7o max 



(t)]}, 



with 



"0(7, T ) = 



7 ~7c 
V 



(22a) 



(22b) 



and 7™ ax is derived below [Eq. (|29a|) ] . 

The other case of interest is m = 1 and arbitrary (d, n) 



C d 



N (-r, T ) = ^TT CX P [(crrVr - $ errV)" +1 j x {9 (7 - 7c) - 9 [7 - ^(t)]} . 



with 



Q'l = erf 



t -t \ J_ 7 - 7c 



(23a) 



(23b) 
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and 7i nax (r) is given by Eq. ([29b)) . 

From the above time-dependent expression we obtain the GSD at full crystallization. For t — > oo, Eqs. 
for m = becomes 



obtained 



AT (7, r — > oo) = 



C d 



Vn 



7 - 7c 



cxp { (-iy 



ti In 



7 - 7c 
V 



n+1 



e (7 - 7c) - e ( 7 - 1) 



The decay of the growth rate is exponential in this case. Of particular interest is the case is when n = d 
The distribution then takes the form 24 



(24) 
1 = 1. 



AT (7, r — > 00) = 



1 



V 



cxp ■ 



7 -7c 

e ( 7 - 7c) - e (7 - 1) 



tfln 



7 ~7c 

Vn 



(25) 



This result is closely related to the lognormal distribution. As shown in Ref. [2J this is seen more clearly by noting 
that, in general, p c <C which implies that for < 7 C < 7 < 1, the equation simplifies to 



AT>(7, r — > 00) 



,io exp {~ $ ln7 ] 2 } 



7 



6(1-7). 



(26) 



The close relation of this result to the lognormal distribution [sec Eq. (|B7[) in Appendix [XT] for the expression written in 
physical units] and its application to full crystallization of amorphous silicon were presented in Ref. [24|. We emphasize 
that obtaining a lognormal- like distribution as a solution of a partial differential equation is a quite remarkable result , 
and we elaborate on its significance in Sec. IIVEI 

Contrary to the case m = 0, no significant simplification of Eqs. (|23[) (m = 1) is obtained in the limit t — > 00. The 
above expressions allow studying the behavior of the theoretical distribution in the following sections, and were used 
to analyze experimental data during,— and at full crystallization 2 ^ of amorphous silicon. 



C. Maximal grain size 

The solution presented in the previous paragraph contains the maximal grain size 7™ ax (T) that can be observed in a 
sample undergoing RNG crystallization. This quantity is obtained from dp — v m (t)dt. Integrating one immediately 
obtains in dimensionlcss variables 

7 r x (r) = ^-M = 7c + - / v m (r')dr' = 7c + u m (r ,r). (27) 



*r J to 



For v m (r) given by Eq. (|13b[) . the above equation can be written in terms of the upper incomplete gamma functions, 
r[a, x], 



7 m "(r) = 1 — ^— r 
7m 1 ' m+1 

Explicit expressions of Eq. (|2"5)) with m = and m = 1 are 



1 / \ m+i 

1 / T - TO 



m+1 



(28) 



_ max 

7o 



7 max (r) = !_ Vo yI erfc IJ! I m=lj (29b) 



(r) = 1 - V a e- (T - To),u , ?7i = 0, (29a) 
■zr erfc (— — 

2 \ t r 

with erfc = 1 — erf. 

It is instructive to digress from the main path of the paper and determine the maximal grain size in the case m = — 1, 
which corresponds to a constant growth rate = evg. We obtain 

7_T x (r) = 7c + eV ( ^) ■ (30) 
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We note that 7" J f x (r — > oo) — > oo. This is consistent with previous work that showed unbounded maximal grain size 
but is not physically justified unless coalescence is taken into account. We do not analyze this case further here, as 
we would have to redefine our dimcnsionless grain size 7. 

One can also determine the maximal grain size p^ once crystallization is completed (t 3> max{i c „, t c i}). In the 
dimensionless formalism used in this section this corresponds to 7^ = 1. The expression for p^ can be transformed 
into an interesting relation between fundamental parameters of the model for any m [see Eq. (|Bllj) in Appendix IB] 



tcvVQ 

This equation can be used as a self-consistency check or to determine the value of one parameter once the others 
have been measured. 



(31) 




FIG. 2: Maximal grain size as a function of time (7™ ax (T) — 7c)/(l — 7c) = (p™ ax (' r ) — Pc)/(Pm — Pc)- The normalization is 
such that the curves do not depend on p„ Vo = t cv vo (see text). The thick lines are for m = and the thin lines m — 1. For 
each set, we have t r = 0.5 (black solid), t r = 0.75 (blue dashed), and t r = 1 (red dotted). The saturation of 7™ ax (r) is shown 
to occur at earlier times with decreasing t r and/or increasing m. 



The time-dependence of the maximal grain size is depicted in Fig. [2] for m = 0, 1 and reasonable values of t r M- Within 
our model w m (r) given by Eq. (|13b[) is proportional to Vo- Hence, if follows from Eqs. (|28[) . (|B10c[) . and (|B11[) that 
(7m ax ( T ) — 7c) /(I — 7c) = (Pm ax ( T ) ~~ Pc) I {Pm — Pc) docs not depend on Vo; the results in Fig. [5] depend only on m 
and t r . It is also worth pointing out that, except for early stages of crystallization, the inequality p™ ax (i) S> p c is 
generally satisfied and Fig. [2] essentially displays the ratio Pm ax ( T )/Pm- The figure demonstrates that the radius of 
the largest grain observed in the sample saturates rapidly in time to reach the size observed at full crystallization 
(calculations of the GSD for the parameters of Fig. [5] show that full crystallization is obtained for At > 8 if m = 
and At > 4 for m = 1.) A decreasing value of t r (e.g., faster decay of the growth rate at constant nucleation rate) 
enhances the rate at which p™ ax reaches p™- A similar conclusion is reached with increasing m. This behavior results 
from the fact that both a decrease of t cv and increase of m accentuate the time decay of the growth rate, which -as 
discussed in previous sections- reflects the impingement caused by neighboring grains. 

The analytical determination of the maximal grain size in terms of fundamental parameters of the model is of 
particular interest for those interested in a probabilistic approach using computer simulation to determine the grain 
size distribution because, as discussed in Ref. [t| the proper account of rare events (the size of the largest and smallest 
grains) is essential for describing the data with the adequate distribution and obtaining the correct value for its average. 



IV. CHARACTERIZATION AND TIME-EVOLUTION OF N(j, r) 

This section is devoted to the characterization of the grain size distribution obtained in Eq. (|21[) . First we analyze the 
influence of the model parameters t r and vq, the ratio of critical times and the growth rate, respectively. This leads 
to scaling properties for vq in the limit r — > 00. We then consider the time-evolution of the distribution. Finally, we 



12 



show how for certain classes of solutions (n, to) the distribution takes the lognormal form in the asymptotic limit of 
large time. 

To proceed with the numerical analysis, it is appropriate to consider the normalized GSD 

^ (7 ' T) = ^P (32a) 

where the total number of grains at time r is given by 

r°° rim^( T ) 
N(t)= / N(*y,r)dj = / N(j,r)d-f. (32b) 

Jo J~f a 

This normalization procedure eliminates the constant factor d defined in Eq. (|15[) . As mentioned earlier, this is the 
only term containing the nuclcation rate la and the results discussed in the remainder of this paper are therefore 
independent of the explicit value of Iq . Such normalization is also useful for comparing the theory with experimental 
dataM 

The figures in the present paper all depict the normalized GSD N(j,t). Thus, the area under the curve is one. It is 
also worth noting for the next sections that none of the distributions discussed here contain divergences. 



A. Moments and maxima of the distribution 



The time-dependence of the distribution can be characterized in terms of its moments. In particular, the first three 
moments, which give the mean p, the variance tr 2 , and the skewness 71 of the distribution, c 2 , 71 are central moments 
and the latter is normalized. The two higher moments give an indication about the spread and asymmetry of the 
distribution about the mean. All moments are calculated for the normalized GSD [Eq. (|32]l ] and thereby independent 
of the nucleation rate coefficient Iq. We define 



<t 2 (t)^ 2 (t) 



7i (r) 



3/2 
M 2 



c?7 7 A(7,r), 
d l [7- K T )} 2 N(-f,T), 

1 f°° 3 _ 

dl [7 - M^)] N(7,t) 



3/2 , 
H J0 



(33a) 
(33b) 
(33c) 



The definitions are given for the dimensionlcss GSD written in terms of 7. Moments for N{p 1 1) are given in Appendix 
151 The conventional notation for the third moment, 71 (always written with its index), should not be confused with 
the variable 7 = (never written with an index.) We calculate the time dependence of the mean and variance 

for the case to = and n = 1,2,3 in Sec. IIVDI 

In some cases, it may be of interest to compare the radius 7^'™ , (i = 1, 2, . . . ) for which the distribution is maximal 
to the mean p, of the GSD. Since it turns out that under certain circumstances the GSD [Eq. (|21|)] has more than one 
maximum (see Sec. IIVD|) . we add the index i to 7™'™. We thus define 7™'™ ^ by 



cW( 7 ,t) 



d-j 



9 2 iV( 7 ,r) 



97 2 



< 0. 



(34) 



The analytical form of 7^'^ m can be determined from the zero of the derivative 



m+l 



For the case n = d = 1 and m = the maximum in the open interval (7 C , 7m ax (i")) is given by 

7^ m= V) = 7c + V [e- 1 ^ + e"(— )/^ 



(35) 



(36) 
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B. Grain size distribution at t — > oo 



For t — > oo, the (unormalized) distribution is given by Eq. (|24|) for to = 0. We remind that in this case the time-decay 
of the growth rate is exponential. The normalized distribution N(-j, r) is shown in Fig. [2ta) for n = d = 1, 2, 3 and 
the parameter values specified in the caption. These are physically reasonable choices of parameters as discussed in 
Refs. [24| and All curves have identical upper cutoff (the upper cutoff is out of the range of the figure for m = 1 ; 



T 




o 

x 



T 



0.2 0.4 0.6 O.i 

y = r/r~ = p/p~ 




y = r/r =plp (x 10 ) 



FIG. 3: Normalized distribution N{p/,T — > oo) for n = 1 (black solid), n — 2 (blue dashed), n = 3 (red dash-dotted line), (a) 
m = 0. (b) m = 1. The abscissa is the normalized grain radius 7 = p/pmt which is independent of n (see text). In this and 
the following figures we consider generic parameters chosen in the range of interest to the experiment^ t r = 0.75, Vo = 1. 



Fig. [3] right). For to = both the mean value and the variance of the distribution decrease with larger n; the peak of 
the distribution is sharper and shifts to lower values of 7 for increasing n. Since we identify n with the dimensionality 
of the system, we expect the GSD of three-dimensional crystallization to be sharper than that of thin films when the 
thickness of the film is smaller than the average grain size at full crystallization. Furthermore, the majority of grains 
have smaller size in three dimensions than in two dimensions. 

The distribution is also affected by the choice of m. Choosing m = 1 (Fig. [2 right) instead of to = (Fig. [2l 
left) implies a faster time-decay of the growth rate. From Fig. [2] we obtain the physically intuitive result that a 
stronger time-decay of the growth rate leads to sharper peaks, with a maximum located at smaller grain radii. This 
is emphasized by the difference in abscissa and ordinate scales between to = and to = 1 in Fig. [3] Replacing 
the exponential decay (m = 0) of the growth rate by a Gaussian decay (to = 1) has a dramatic effect on the grain 
size distribution. This results in the unambiguous choice to = to describe the experimental data of solid-phase 
crystallization of amorphous siliconi 24 ' 33 



C. Dependence of JV(7,t) on t r and Vo 

The grain size distribution -/V(7, r) is a function of d, n, to, t r , to and vq. However, not all parameters need always be 
known to determine the GSD. Considering Eqs. (fT8|) . (|28|). and (|B13[) . we note that it is possible to reduce the number 
of parameters that appear in Eq. (|2ip . For example, from Eqs. (|18p and (|28p . we can express <7 m (7, t) in terms of t r , 
7c 7™ ax (7"), and To- Thus, the normalized GSD can be written in terms of the latter parameters, and Vo and lo are not 
explicitly needed in the expression. Conversely, it is possible to express the GSD in terms of the two latter quantities, 
thereby, removing other parameters. For example, for to = 0, Eq. (|B13|) implies Vo = (1 — 7 c )r _1 [(to + 2)/(m + 1)], 
which can be used to write (7 - 7 C )/V = T [(to + 2)/(to + 1)] (7 - 7 C )/(1 - j c ) in Eqs. and ([2"3|). As is often the 
case, 7 C <C 1, which implies that the ratio is essentially 7. Then, in the limit r — > 00, the GSD [Eq. (|24p] only depends 
on t r and Cd and the normalized GSD only on t r . At finite times 7™ ax (T) remains present and is directly proportional 
to Vo- However, 7™ ax can be determined experimentally, and Vo is not required. To summarize, the choice of which 
parameters are needed and which can be obtained from the expressions above depends on the particular situation 
under consideration. A natural choice of parameters to discuss the properties of the normalized distribution is t r and 
Vo- We discuss in this section the general dependence of N(j, r) on these two parameters. 
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To study the generic influence of t r on the behavior of N(-f,r), it turns out to be sufficient to consider Eq. (f2Tj) 
in the limit t — > oo when full crystallization is achieved. This limit has been calculated analytically for classes of 
solutions (d,m,n) relevant for experimental studies in Eqs. (|22M24j) . Figures 0] and [5] highlight the influence of t r and 
vo, respectively, on the shape of the distribution for m = (left column) and m = 1 [right column; the rows are for 
n = 1 (top), n = 2 and n = 3 (bottom)]. The first general observation is that in all cases, the distribution displays 




FIG. 4: iV(7,r — > oo) for different critical time ratios t T = \/t cv /t c i = 0.5 (black solid line), t r = 0.75 (dashed blue line), and 
t r = 1 (dash-dotted red line). In Figs. (a,b,c) the magenta dotted line is for t r = 0.6. In Figs. (d,e,f) the magenta dotted line 
is for t r — 0.89 and the green dash-double dotted line is for t r = 0.92. The left column is for m = and the right column for 
m — 1. Rows are for n = 1 (Figs. a,d), n — 2 (b,e), and n = 3 (c,f). Note the different scales on the abscissa and the ordinates 
between the two columns. For all figures Vo = 1- 



one maximum, and cutoffs at 7 — Pc/Pm an d 7 = 1. While the latter is obvious from Eq. (1211) . the former is only 
true at large times within or model, as will be seen in the next section. 

As t r = t cv /t c j < f decreases, the number of small grains increases at the expense of the formation of larger grains. 
This is confirmed in Fig. 21 which shows that irrespective of the value of n and m a decrease of t r results in an 
increase in amplitude, a sharpening of the peak, and a shift of the maximum to smaller values of 7. Furthermore, 
the properties observed on Fig. [3] are also found in Fig. |4] Thus, the qualitative features inferred from Fig. [3] do not 
strongly depend on the particular value of t r when taken within a physically reasonable range. It is interesting to 
observe that the position of the maximum of the distribution for n = 2, 3 and m = 1 is essentially insensitive to the 
value of t r unless the latter is very close to one. Even then, comparing the order of magnitude of abscissa and ordinate 
scales of Figs.EJe) and|4jf) and|4ja) and|U[b), the peak of the distribution barely shifts for m — 1 and increasing t r . 
On the other hand, the amplitude of the maximum strongly varies with t r . Finally, the distribution for n = m = 1 
and t r = 1 is rectangular [lower red dotted line in Fig. Hfd)]. This result is derived in the next section. 
We now study the influence of the growth rate amplitude Vo on the grain size distribution. We can also limit our 
considerations to the limit t — > 00. The influence of the growth rate amplitude Vo is shown in Fig. [5j Physically, an 
increase of Vo implies a faster growth of the grains in a given time interval, thus, leading overall to the presence of 
larger crystallized grains in the sample. According to Fig. [4j the ratio t r determines how sharp the peak is at given 
Vo- At fixed value of t r and for all n, an increase of Vo results in a decrease and broadening of the peak as well as 
a stretching of the distribution to higher values of 7. As in Fig. HI the peak for m — f is also barely shifting with 
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FIG. 5: JV(7, t — > oo) for different growth rate amplitudes Vo- Black solid line: Vo = 0.5 (left column) and 0.6 (right column). 
Dashed blue line: Vo = 1. Dash-dotted red line: Vo = 1.5. Left column: m = 0. Right column: m — 1. Figs. (a,d) are for 
n — 1, (b,e) for n = 2, and (c,f) for n = 3. In Figs. (a,b,c) t T = 0.75, and in Figs. (d,e,f) t r = 0.9. 



change of Vo, as opposed to the case m = 0. Although a trend similar to Fig. [4] is observed when increasing n, the 
growth and shift of the peak are less pronounced. Note that the blue dashed curves with parameters t r = 0.75 and 
Vo = 1 are the same in Figs. @] and [5j which allows a scaled comparison between of the two figures. 
To conclude this section, Figs. 0] and [5] show that the behavior for m = 1 as a function of t r and Vo is similar to that 
for m = 0, but the peaks are much sharper, larger in magnitude (we remind that none of the distributions contain 
divergences), and located at smaller values of 7. In addition, we note that the shift of the peak with increasing value of 
t r or Vo is much less pronounced for m = 1 than for m = 0. This reflects the difference between the exponential- type 
decay (m = 0) and the Gaussian decay (m = 1) of the growth rate. It also reflects the presence of the logarithm or the 
error function of 7, respectively, in the distribution. The different behavior of the distribution for m = and m = 1 is 
distinctive enough to determine which distribution is most suited for a set of experimental data. For example, in our 
analysis of crystallization in amorphous Si thin films (Refs. and |33[). one is unambiguously led to choose m = 0. 
Remember finally that in the KAMJ model, n = d is the dimensionality of the crystallization process. Thus, the fit of 
experimental data for different values of n could be used to determine the dimensionality of the crystallization process. 



D. Time-evolution of the distribution 

In the previous section, we discussed how model parameters influence the grain size distribution in the fully crystallized 
limit t — > 00. The results qualitatively hold for all times To < r < 00 as well. We now study how the distribution 
evolves in time for fixed values of the parameters as the RNG process takes place, starting from an amorphous solid 
in d dimensions. 

It is instructive to consider first two simple cases of the general solution. When the critical times for nucleation and 
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growth are identical, t cv = t c j = t c (implying t r — 1), Eq. (f2T|) reduces to 

iV( 7 ,r) = -^e^-")^ 1 -^- o)" +1 ] x {e( 7 -7c)-e[ 7 -7^(r)]}. (37) 
Assuming further that n = m, we obtain the simple form 

N fr> T ^^TT X - 7c) - e [ 7 - 7 r x (r)] }■ (38) 

This latter case is depicted in Fig. IDJd) (red dash-dotted bottom line, close to the abscissa) and approximates for 
arbitrary values of n and m the grain size distribution at times very close to the incubation time io-~ However, our 
description of the nuclcation process [Eq. (QJ] is rather coarse and is not expected to reproduce accurately the early 
stages of crystallization^ 

Comparing the above two expressions, we note that in the first case the distribution displays an exponential dependence 
of time [remember that a m — cr m (7, r)], whereas in the second case the distribution is constant in time. These two 
examples emphasize the fact that the exact functional form of nucleation and growth rates, and, in particular, their 
time dependence [the choice of n and m in the power laws of Eq. (|13[) ]. is critical in determining the grain size 
distribution N(p,t) at all times, including in the asymptotic regime t — > oo. 

To analyze the time dependence, we focus on the classes of solution defined by m — 0, 1 and n — 1,2, 3, where analytical 
solutions of Eq. $Q can be derived, and that also appear to be most relevant for the description of experimental 
data. 24 ' 33 The expressions for the grain size distribution are given by Eqs. (|2"2"j) for m — and Eq. (f2"3"]) for m = 1. 
Figure [6] displays the time-evolution of the distribution in three-dimensional plots for the various values of m and n. 




FIG. 6: Time evolution of the distribution JV(7, r) for n = 1 (top left), n — 2 (top right), and n = 3 (bottom). m = 0. t r = 0.75, 
Vo = 1. The lines at constant At = t — to are depicted in Fig. [7] left column. The curved line in the (7, At) plane shows the 
maximal grain size as a function of time and corresponds to the curves of Fig. [2] 



Three observations can be made on Fig. [6] First, considering the intermediate to asymptotic time development when 
growth dominates over nucleation an increase in the value of n leads to a more pronounced peak and decreasing value 
of the mean. Thus, in accordance with the description of Figs. [3] and [5l the faster the decay of the nucleation rate 
the smaller are the grains in average once crystallization is complete. Second, the specific time-dependence of the 
nucleation and growth rates has a notable influence on the early stage of crystallization. Whereas for n = 1 (slow 
decay of the nucleation rate) the distribution has one broad peak at low 7 that further broadens and decreases in 
magnitude over time, for n = 2 a sharp peak is observed at small radii of the grains (henceforth called the "nucleation 
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FIG. 7: N(fy,r) at different times Ar st-to (cuts from Fig. [6]). Left column: m = 0, right column: m = 1. Top row is 
?i = 1 for (a) and (d), n — 2 for (b) and (e), and n — 3 for (c) and (f). Note the abscissa and ordinate scales, especially for 
Figs, (e) and (f). The time dependence is discussed in the text. t r — 0.75, Vo = 1. The curves displayed for m = (left 
column) correspond to the thin lines at constant Ar in Fig. [6] 



peak") and a second smoother peak at larger values of 7 (called the "growth peak"). The case n — 3 is similar, but the 
two peaks are better resolved. This is discussed in Fig. [7J Figure |5] also displays a line delimiting the distribution at 
large 7. This line corresponds to 7™ ax (r) and was depicted in Fig. [2] Note that the line saturates as r — > 00 because 
of the time dependent growth rate. This contrasts with the result found in previous work, where 7 max (r — > 00) — > 00. 
Figure [7] displays the cuts shown in Fig. [6] for specific times r. The case (n, m) = (1,0) (exponential and Gaussian 
decay of the growth and nucleation rates, respectively), Figure UJa), has the shape of a lognormal distribution at 
infinite time. In this case, the grain size distribution has little structure at early stages. The case m = 1 = n 
[Eq. ([37])] is shown in Fig. [TJd) and corresponds essentially to the product of two Gaussian functions with different 
prcfactors that lead to the competition between a Gaussian time-decay of nucleation and a Gaussian saturation of 
the grain growth. Figure [T^d) displays the tail of the Gaussians. 

More interesting are the cases n = 2,3. In Figs. [T^b) andJTJc) (to = 0) the sharp nucleation peak near j c (typically 
7 C ~ 10 -3 — 10~ 4 ) and the broad growth peak described in the three-dimensional plots are clearly resolved. The time- 
evolution of these two peaks can be interpreted as follows. At early stages nucleation dominates the crystallization 
process leading to the formation of a large number of nuclei in a short time, and thus to a sharp peak at small radii. 
As time passes the nucleation peak is depleted and extends to larger radii [case Ar = 1 on Fig. [7fb)]. Because the 
decay of nucleation is super Gaussian an intermcdiatc-radii peak emerges as a result of the slower growth decay. The 
increase in the growth peak at intermediate values of 7 occurs at the expense of the nucleation peak (near j c ). The 
front of the distribution shifts substantially as the broad peak starts emerging, then moves a little as the nucleation 
peak decreases, and finally shifts again as only growth determines the distribution. 

Figure [5] displays the time dependence of the mean and the variance for the case to = and n = 1, 2, 3 (left column 
in Fig. [7]). We note that for n = 1 both the mean and the variance increase smoothly to reach the value at infinite 
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FIG. 8: Time-dependence of the mean fi(r), Eq. ()33ap , and 
figures are for m — 0, and n = 1 (solid black), 2 (blue dashed) 
for the time-evolution of the grain size distribution depicted ii 




AT= T- T q 

variance er 2 (r), Eq. (|33b|) . of the grain size distribution. Both 
, and 3 (red dash-dotted line). The mean and variance are those 
l Fig. left column. t r = 0.75, Vo = 1. 



time. On the other hand, the cases n = 2 and 3 display a "kinklike" structure in the range At <S [1,3]. This structure 
reflects the fact that the nucleation peak decreases and broadens at first but remains very sharp, until At « 1. The 
appearance of the growth peak leads to an abrupt and large increase in the mean and the variance, which proceeds 
until the nucleation peak has completely disappeared. Once this stage is reached, both the mean and the variance 
increase marginally to their values at full crystallization. 

The experimental observation of two peaks in the GSD is difficult and may often not be possible. For example, in 
the solid-phase crystallization of amorphous silicon only one peak is resolvedi 24 ' 33 There are two possible reasons for 
this. One is that the nucleation peak is rather sharp and may be below experimental accuracy. The other is that the 
two peaks strongly depend on the time-evolution of nucleation and growth rates, but also on the early stage dynamics 
of nucleation. In the present model, we described the source term by a Dirac delta distribution [Eq. (j¥])], which is 
a rather crude approximation of real systems. Instead, the delta distribution should be replaced, for example, by a 
Gaussian of finite variance, which spreads the formation of nuclei over 7^ This modification broadens the nucleation 
peak and is expected to fill the dip between the nucleation and growth peaks and may prevent the resolution of the 
two peaks. 

It is interesting to note that a Gaussian depletion of the growth rate obtained for m = 1 also leads to the disappearance 
of two well- resolved peaks. This is depicted in Figs. [7(e) and[JJf). Nonetheless, the dynamics of crystallization is 
otherwise very different for m = and m = 1. One essentially observes a tall peak at small values of 7, the magnitude 
of which decreases and broadens in time to become the infinite time distribution. But the maximum of the peak for 
m = 1 is orders of magnitude larger than for m = [note the different abscissa and ordinate scale of (e) and (f)]. In 
addition, because of the more rapid decay of the growth rate, the second peak has its maximum at smaller radius and 
is consequently fully absorbed by the nucleation peak. This is also the reason why the peak maximum barely shifts 
with time. 

Figures [Jj and [8] convey the message that the time dependence of nucleation and growth rates is essential for the time- 
evolution of crystallization processes, and also leads to very different shapes of the final GSD at full crystallization. 



E. Discussion of the general solution and the derivation of a lognormal distribution 

We discuss the results of the previous paragraphs in the context of work done on random nucleation and growth 
processes and more generally work involving lognormal distributions. We showed in Ref. [24| and Eqs. (|2"5j) . (|2l)]) . and 
(|B7[) that a lognormal distribution is found at infinite time (more precisely, when t 3> t c i, t cv .) This distribution 
has been so widely studied for many years that it may seem surprising that it has not been previously derived from 
a partial differential equation. To the best of our knowledge, this is, however, the case and is certainly true in the 
context of random nucleation and growth processes. Although observed in many instances, there are only few recorded 
attempts to actually derive this distribution from first principles. This is not to say that the lognormal distribution 
has never been derived. 

In fact, studies involving the lognormal distribution can be divided in two categories. One category, which in- 
cludes the majority of publications on the subject, postulates a lognormal distribution (or more involved varieties 
thereof) ^ ] 30 ] 32 i 36 i 40 A fit to experimental data is then performed to determine the parameters of the distribution. 
Work belonging to this category does not derive the lognormal distribution and can, for example, not express the 
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mean and variance in terms of fundamental parameters of the system. 

The second category forms a much smaller group of publications, in which a derivation of the lognormal distribution 
is proposed. In an insufficiently cited paper,— A.N. Kolmogorov presented what is likely the first derivation of the 
lognormal distribution, applying probabilistic and statistical arguments to crushed powders, where previous experi- 
mental studies had shown that the lognormal distribution delivers a good fit to the data. Shortly thereafter, Epstein 
discussed the same grinding process and presented an alternative, somewhat clearer derivation,— demonstrating that 
under fair ly g eneral conditions the lognormal distribution is obtained as a result of the central limit theorem (see 
also Refs. Eol and l9h . It is nowadays widely spread in the literature that a lognormal distribution can be obtained 
for multiplicative random processes (as opposed to additive processes, such as the random walk) using the central 
limit theorem applied to the natural logarithm of the product of probabilities. The limitations and problems of this 
approach have been discussed in detail in Ref. [|1 together with their remedies. 

Our derivation differs from the work of Kolmogorov and Eptsein in several ways. First, our study of RNG processes 
starts at the opposite end of particle formation. Instead of breaking larger particles into smaller pieces, RNG starts 
with the creation of nuclei that grow over time. Thus, our description of the grain size distribution must include two 
physical phenomena: a source term accounting for the initial formation of grains and the growth of these grains. 
References [25| and [26| do not contain any information on the source term (except for a postulated initial distribution) 
and how it affects the time-evolution of the particle size distribution. Second, we derive the lognormal distribution 
from a partial differential equation rather than probabilistic arguments. This is thus done in the same spirit as 
Avrami's approach to the kinetics of phase transformation, although the latter only determined the fraction of 
transformed material and not the grain size distribution discussed here. We do not refer to the central limit theorem 
in the derivation and overcome thereby limitations related to this approach^ On the other hand, we explicitly use 
the conservation of mass, which is implicit in Kolmogorov and Epstein's derivations. Finally, contrary to the latter 
authors, our derivation enables us to write an explicit closed analytical expression for the grain size distribution at all 
times, and not only in the asymptotic limit. This allows us to study the time-evolution of the grain size distribution 
observed in the crystallization process of amorphous Silicon (see Refs. l24l andl33l) and other materials. We also stress 
the fact that our derivation is done for arbitrary dimension of the grain formation and the truncated lognormal 
distribution (together with the caveats mentioned in Sec. Ill Bl and Appendix |B|) is only one of the solutions obtained, 
namely in the case n = to + 1 = 1 and in the asymptotic limit of large times. 

Another remark concerns the partial differential equation and the determination of its analytical solution. Partial 
differential equations similar to the one established in Ref. for RNG processes and discussed here have been 
studied at length in the literature^ They have been applied to describe various phenomena in nature and, in some 
cases, an analytical solution has been derivedi 11 ' 12 ' 35 ' 36 However, none of the papers have recognized and proven that 
the lognormal distribution is actually a solution of such a PDE. Our derivation clarifies why this has been overseen. 
Although a formal solution of the PDE can be obtained by quadrature, specific additional conditions on the functional 
form of v(r, t) and 2?(r, t) are required to derive the lognormal solution. Generally, when specified in the literature, 
polynomial forms in r and simple exponentials in time have been considered for these functionals, and these do not 
lead to a lognormal-type distribution! 11 ' 35 i 36 The class of PDE considered here involves more complicated exponential 
forms of the time dependence that result from the physical analysis of RNG processes^— and these do lead to the 
obtained lognormal-type forms. In this context, it is interesting to note that the closest expression to a lognormal 
distribution is obtained for d = 1, where the functional form of ~v(r,t) is exponential in time (to = and v has no 
dependence on r), whereas T>(r,i) is Gaussian in time (n = 1). This leads to the conclusion that when a lognormal 
distribution is experimentally observed for grains described in terms of a scalar such as the volume or the average 
radius of the grains it is indicative not only of a particular size dependence of the source-and-growth terms but also 
of the dynamics of crystallization. The general conditions under which a lognormal distribution is a solution of the 
PDE is a topic of its own and will not be discussed further here. 

Finally, the results obtained in this paper are not only interesting in themselves, but allow also to go one step further 
in the phenomenological description of RNG processes that may be useful for applications, since we determine all 
features of the distribution in terms of the fundamental parameters of the RNG model, namely, Iq, Vq, t cv /t c i, to, 
and p c . The theory may, in principle, even be applicable to other phenomena not limited to physics and where 
distributions display, for example, a lognormal-like behavior asymptotically. As mentioned earlier, it is well known 
that such distributions are observed in biology, economics, sociology, etci 30 ' 32 Clearly, not all these phenomena can 
necessarily be described in terms of the differential equation But, it cannot be excluded that for certain problems, 
where a conservation law is present, and v m (r, t) and T>(r,t) have appropriate context-given meaning, the equation 
and solutions discussed here may be relevant. It is beyond the scope of this paper to discuss such generalization. 
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V. CONCLUSIONS 

The paper presents a description of the time-evolution of the GSD in the crystallization of a dimensional solid. 
We provide an analytical derivation of the general solution of the partial differential equation for RNG and the 
Kolmogorov-Avrami-Mehl- Johson model for the fraction of material available for further crystallization. The general 
analytical expression for the GSD, Eq. ((21]) with Eqs. (fT"6|) and (|10p (see also Appendix IB1 for expressions that can be 
used for experimental data analysis), has been divided into classes defined by (d, n, m), where d is the dimensionality of 
the growth process, n and m determine the time-decay rates of nucleation and growth, respectively. A key ingredient 
of the present theory is the introduction of an effective time-dependent growth rate in addition to the usual effective 
time-dependent nucleation rate. This new time dependence takes into account the fact that impingement inhibits the 
growth of grains, which reduces the actual average growth rate in time. 

The main conclusions can be summarized as follows. First, the time decay of nucleation and growth rates plays a 
major role for the evolution of the distribution during crystallization. The time decay strongly affects the shape of 
the distribution, even at full crystallization. Thus, the time path followed by the system to reach the equilibrium 
distribution is essential. 

Second, in the particular case of an exponential decay in time of the effective growth rate (m = 0), the grain size 
distribution develops asymptotically in time to a form that is very close (n = 2, 3) or even essentially identical (for 
n = 1) to the well-known lognormal distribution. The only difference being the presence of cutoffs at small and large 
grain sizes and an overall constant prefactor [see Eqs. (|26|) and (|B7j) ]. This result was pointed out in Ref. Hfor d = 3. 
Its extension to any dimension, and its significance are discussed in the present paper. It is remarkable that the 
derived GSD has a shape that is closely related to the lognormal distribution. It would be of interest to define more 
precisely the conditions on the functional form of source-and-growth terms under which a lognormal-likc distribution 
is a solution of the PDE. Note that our study does not include possible additional contributions to the crystallization 
process. For example, coarsening is explicitly discarded from our considerations, although this phenomenon is also 
known to lead to a lognormal grain size distribution in the asymptotic time regime. This, however, is known for a 
long time. The realization that coarsening is not necessary, as discussed here and in Refs. y@ and|H, is a more recent 
insight on crystallization. 

Third, the relative time decay of nucleation and growth rates affects drastically the mean and the variance of the 
distribution. For the model used here to describe nucleation [a Dirac delta distribution, Eq. Q], the model even leads 
to the presence of two peaks at intermediate stages of crystallization. As discussed in Sec. IIVDI the dip separating 
the two peaks is likely to disappear when replacing the Dirac distribution by a Gaussian nucleation rate with finite 
variance. This requires further work. 

Finally, the presence of a time dependent effective growth rate leads to a finite maximal grain size at all times, 
including for t —> oo [Eqs. (|28p or Eq. (|B10c[) ]. This contrasts with existing theories that have a constant effective 
growth rate and for which the distribution does not have a maximal grain size above which the distribution vanishes. 
One interesting prediction of the model is the dependence of the distribution on the dimensionality of the crystallization 
process. For example, it is of interest to study the grain size distribution of thin films with various thicknesses. For 
film thicknesses smaller than the typical grain size, the experimental distribution should be best described with d = 2. 
As the thickness of the film increases beyond the average grain size, the theoretical expression with d = 3 should offer 
a better description of the data. Also, our theory provides a fundamental relation between parameters of the model 
and experimental quantities [Eq. (|3ip]. 

An example of application of the theory to the time-evolution of the GSD in solid-phase crystallization of amorphous 
silicon has been presented in Refs. [24| and l33l . We emphasize that contrary to many studies that rely on an ad hoc 
lognormal distribution to describe experimental data, we propose a description based on the creation and growth of 
grains. Thus, the shape of the distribution is derived and not postulated and relies on the knowledge of a few physical 
parameters describing the system. The derivation of the GSD has a predictive power for the physical description of 
the time-evolutio n 24 ' 33 that empirical approaches cannot offer. 

Finally, the general formulation proposed in this paper for processes involving a source-and-sink term and a growth 
term should be applicable to a wider range of phenomena, where the grains are other entities, and where source- 
and-growth terms need to be defined appropriately. A few examples were offered in the introduction. We hope that 
the formulation and the derivation of the solution are written in general enough terms to initiate its application in 
other fields. We emphasize, however, that not all crystallization phenomena can be described by the KAMJ model. 
Therefore, we do not expect the solutions derived in this paper to be equally applicable to all materials. When 
they do, it is an indication that the Kolmogorov-Avrami-Mehl- Johnson model for the time-dependent fraction of 
the transformed phase is a good description of the crystallization process, as can be experimentally observed for the 
crystallization of amorphous silicon^ When they do not, the original differential equation should be considered again, 
with other forms of the nucleation and growth rates, possibly including anisotropy and size dependence. 
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Appendix A: Table of variables, parameters and relations 



We provide a table (Table [TJ summarizing all variables, parameters and their definitions. In addition, we present 
below the table some relations between parameters. 



TABLE I: Definitions of variables, constants and functions used in the 
paper. 

N(r,t), N(t,t) Grain size distribution (GSD), Eq. ©, and normalized (GSD) defined in Eq. ([32)) . 

N(r, i) Auxiliary dimensionless function introduced in Eq. (|12[) to transform the partial dif- 

ferential equation (PDE) in a form that can be solved analytically. 

d Dimension of the crystallization process. Although the calculations are valid for any d 

we consider specifically d — 1, 2, 3. 

r = (r 1 , . . . ,rd) d— dimensional vector of magnitude r, the components of which are the semi-axes 

n > T2 > ■ ■ • > rd of the ellipsoid that models a particular grain. 

p Radius of the spherical grain. The latter is the limiting case of an ellipsoid with all semi- 

axes rj — p (j — 1, . . . , d) and the previous d— dimensional vector reads r = (p, . . . , p). 
Consequently, r = \fdp. 

p c Critical radius. Radius of the nucleus, the smallest grain that can be found in the 

system. 

p™ ax (t) Radius of the largest grain found in the sample at time t. The radius of the largest 

grain depends on the growth rate. Accordingly, the maximal grain size has an index 
m, which speficies which growth rate is considered, Eq. (|10bl) . 

Pm Radius of the largest grain found in the sample at t — > oo (full crystallization). 

7 = r/r^ — p/pm Dimensionless variable for the size of the grain. Except for Appendix[B]all calculations 

are done using this dimensionless variable. In the text we omit the index m for clarity. 

t, to Time variable t and incubation time to. 

tcv, td, t r critical time of decay for the effective growth and nucleation rate, v(t) and I(t) in 

Eq. (nnj). t r = s/tcv/td- 

t = t/i/t cv tci Dimensionless time variable, to is the corresponding dimensionless incubation time. 

n , m Index of nucleation and growth laws. Eq. (|10|l . respectively. In the paper we focus on 

n = d and m = 0, 1. 

Jo, vo Constant microscopic nucleation and growth rates. 

To , Vo dimensionless constants defined in Eq. ()llb|) . 

Y n (t ) , (t ) , YZ (t ) KAMJ fraction of available space for crystallization and functions determining the time 

dependence of the nucleation and growth rates I(t) and v(t), respectively. See Eqs. ([9]) 
and (flQ) . 

v(t) Time dependent nucleation and growth rates, Eq. (|10[) . 

SI Volume of an ellipsoidal grain defined in Eq. ([5]) 

A Ct d, Aoo^d Surface of the nucleus and the largest grain at full crystallization, defined near Eq. ([7]) 

and l|12|l . respectively. 

Cd Constant prefactor of the grain size distribution, defined in Eq. (|15[). 

o"m,i(7i r ) Function appearing in Eq. (|14l) and defined through Eq. (|16l) . 

p(r) , u 2 (t), 71 (t) mean, variance and skewness of the grain size distribution. They are defined in Eq. (|33p . 

Note that 71 cannot be confused with the dimensionless size variable 7 which is never 
written with the index m (see definition above.) 



A word of caution is necessary with respect to the definition of the rates Iq and vq. In expressions where quantities with 
dimensions (as in Appendix [B]) are used, the definition of Iq an d depends on which variable is used to characterize 
the size of the grains r, p, or the diameter g = 2p, as in Ref. |24| . For example, the growth rate v(t) appearing in 
terms of the different variables is related by if = dr/dt = Vddp/dt = \fdv p and v(t) = d/y/dr is given by Eq. CE3b}. 
Throughout the paper we omit the index r, p or g for the rates as the expression to use is univocally determined by 
the context. A similar argument applied to Jo gives Iq = Jq. It turns out that with this precaution in mind, the 
expressions derived in the paper are essentially invariant with respect to the choice of variable, except the terms in 
the Heavisidc function. This remark about the choice of variables is only relevant for Appendix [B] below since the 
rest of the paper is written in dimensionless quantities and involves only t r , Vo, and X only. 
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Appendix B: Main results written in quantities with physical units 



The distributions derived in the main part of the paper were presented and analyzed using dimcnsionless quantities, 
7 = r l r m = P/Pm f° r tric g ram size and r = t/y/t cv t cI for time. Since one of the goals of the paper is to provide 
expressions for the grain size distribution that can be applied to the analysis of experimental data, we provide in this 
section the main results in quantities with physical units, using p for the radius of the spherical grain and t for the 
time. The expressions below take the exact same form when replacing the radius p by the grain diameter g, under 
the condition that appropriate definition of vq and Iq are considered (see discussion in the previous appendix) . 
According to Eq. (TH|) or Eq. (f2~Tj) , N(p,t) is defined in units of number of grains per unit volume, per unit length (of 
the radius of the spherical grain) or, equivalently, TV dp has the units number of grains per unit volume. Thus, Iq is 
expressed in number of grains per unit length and unit time while vq is in length per unit time. 
For arbitrary n but m = 0, Eqs. (|2"2")) read as 
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The latter is the maximal grain size found at time t during crystallization and at t — > oo for m = and arbitrary n. 
Note that one often has pg° 3> p c - 

For infinite time (full crystallization), the above expression simplifies to 
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Choosing n = d = 1 (m = 0) and t — > oo, one obtains the expression derived in Ref. 
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As discussed in the main text, this distribution is of the lognormal type as seen most conveniently by considering the 

case p c <C p < /0§° 
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where /iog(p) is the lognormal distribution with M as the mean and s as the standard deviation of the variable's 
logarithm. Note that the range of finite values of the distribution [ Eq. (|B7[) ] is (0,pg°), while for the standard 
lognormal distribution it is [0, oo). This physically justified limited range of finite values, together with the presence of 
the multiplicative constant C\ og in front of the expression lead to term the derived distribution as being lognormal-likc. 

Consider now the case m = 1. For arbitrary n [Eqs. (|23|) ]. we have 
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At infinite time, this expression reduces to a form that cannot be readily cast into a lognormal-typc distribution, as 
discussed in the paper. 

Finally, the general expression, Eq. (|21j) for any n and m reads as 
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p = p c + u m (a m ,t) = p c + v m (t')dt' 



where tr m (p, t) is solution of 



with v m (t) defined in Eq. (|10b[) . and p™ ax (£) is obtained in a similar way as Eq. 
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The maximal grain size p^ = lim r ^.oo PnT x ( T ) obtained once crystallization is completed (t m&x{t cv , t c j}) imme- 
diately follows from the above expression 



Hm p2 ax (T) = p% = p c + p™v r 



m + 1 



= Pc + t cv Vo r 



m + 1 



For example, 
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pT 



p c + t cv v a , m = 0, 



= p c + -r- t cv v , m — l. 
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This leads to an interesting relation between fundamental quantities of the model. Reformulating Eq. (jBlip . one 
obtains Eq. (|31[) for any non-negative integer m. This relation quantity can, for example, be used to determine the 
critical time t cv once p^ and v have been measured, and taking into account the fact that, in general, the inequality 
Pm ^ Pc holds. One can also write Eq. (|B11[) in dimcnsionlcss quantities 



i = 7c + v r 



m+1 



(B13) 
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This has been used to obtain Eq. (|28|) . 

To compare the above expressions with experimental data, it is best to normalize the distribution 

*<A*) = l^p ( B14a ) 
where the total number of grains at time t is given by 

r°° rCW 
N(t)= / N{p,t)dp= / N{p,t)dp. (B14b) 

JO J Pc 

This removes the coefficient 

In Sec. IIV Al we defined the three first moments in terms of the dimensionlcss radius 7 = p/p^ for the normalized 
GSD. These quantities are correspondingly defined for GSD in terms of p and we only write the relations between the 
two definitions: 

M*) = /O, (B15a) 
<r 2 p (t) = (M2, P (t) = 0C) 2 M2, (B15b) 
7i,p(i) = 7i- (B15c) 

The definitions in terms of p are supplied by the index. Note that if one considers unnormalized quantities, these 
equations have to be modified accordingly. 

Appendix C: Solution of Eq. © and proof of Eqs. (f20l) 
1. Solution of Eq. © 

We solve the partial differential equation ([7]) and derive Eqs. (fH)) and (|2~Tj) . Using the Laplace transform, 

N(k,T) = C N(j,t) = / AT (7, r)e- fc7 d 7 , (CI) 



with fceR, Eq. ([7]) transforms into a first-order ordinary differential equation in time 

■^-N(k, t) + A « m (r) N(k,T)=t r I n (r)e- k ^ . (C2) 

where v m (r) and I n {T) are given by Eq. (|13[) . The solution of this equation can be immediately found to have the 
form 

N(k, t) = t r e- k ~>" ( I n (a)e- ku ^ a ' T Ua, (C3) 

J T 

where we used the definition ([T7| . The inverse Laplace transform of Eq. (|C3|) leads to the following form of A (7,7-): 

N^,T)=t r l f da e~ A '^ 5(j-B m (a,r)), (C4a) 

with 

= t n r +1 (a - T ) n+1 , B m (a, t) = lc + u m (a, r) (C4b) 

This result was obtained by exchanging the order of the integrals over k and a. This operation imposes convergence 
conditions on the integrand, which arc satisfied for the KAMJ model for nucleation and growth. 
We now use the relation 

S (7 - B m (a, r))=J2 i o t ia ~ "7 } , (C5a) 

l I da V > / I CT =cr m ,■ 
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where we sum over all solutions u m ,i of the equation 

7 = B m (a mti ,T), 



(C5b) 



with J^B m (a, t)| ct 7^ 0. This is Eq. (|16p . Note that cr„ v ; = <7 TOi j( 7 ,T). For the KAMJ model, the sum over i 
reduces to one term because the function B m is monotonous in a for arbitrary values of r and m. Furthermore, the 
condition on the derivative of B m is satisfied for a m .i < 00, which is always true. 

Taking into account the fact that Eqs. (|C4[) are finite when a m ,i G [to,t], we obtain the general solution of Eq. ([7]) in 
the form 



A( 7 ,r) = t r l J2 



■6(r- <T m ,i) <d{a Tl 



To), 



which for the KAMJ model becomes 



t n+L (<T 



ro) 



n+1 



@(t - (J m .i) Q(cr r , 



ro) 



(C6) 



(C7) 



Equation (|14j) follows from this expression. The expressions for a mt i(j,T) are determined from the solution of 
Eq. (|C5b[) and explicit expressions for m = 0, 1 were determined for the KAMJ model in Eqs. (fn 



2. Proof of Eqs. (J20I 

To prove that Eq. (f2~Tj) is equivalent to Eq. ([14")) , we have to prove Eq. (|20p . The product of Heaviside functions on 
the left-hand side of Eq. (f20| is equivalent to the statement r > ct„ v ; > to > 0. Using Eq. (fl~7|) and noting that if the 
expression for v m is non-negative [as is the case of Eq. (| 13b>|> ] . the previous condition on a m ,i implies u m {p m ^,T) > 0. 
From Eq. (|16p . we immediately conclude that 7 > 7c , which is known to be true by definition since 7 = p/Pm- 
On the other hand, for non-negative v m the inequality r > <r m ^ > To > also implies < u m (o m ^,T) < u m (ro,r). 
Using again Eq^. (fT6)) . this condition can be written as 7— j c < w m (ro, r) or, with Eq. (|27|) . 7 < 7c+w m (ro, r) = 7m ax (T). 
This proves Eq. (|20p and, consequently, with Eq. (|C7p . completes the derivation of Eq. ([2~T|) . 
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